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Abstract 

The constraints on an R-parity violating supersymmetric theory arising from the recent 
SuperKamiokande results on atmospheric neutrinos are studied, with special reference to a 
scenario with bilinear R-parity violating terms. Considering both the fermionic and scalar 
sectors, we find that a large area of the parameter space is allowed, in terms of both the lepton- 
number violating entries in the superpotential and the soft R-violating terms in the scalar 
potential, and that no fine-tuning is required. However, the need to avoid flavour changing 
neutral currents puts additional restrictions on the theory, requiring either the R-violating terms 
in the superpotential to be smaller than the R-conserving ones, or a hierarchy in the R-violating 
parameters for different lepton flavours in the superpotential. 
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1 Introduction 



Ever since the evidence in favour of neutrino oscillation has received reinforcement from the Su- 
perKamiokande(SK) results on the atmospheric f^-defhciency Q, intensive discussions have taken 
place on 'new physics' that can give rise to neutrino masses and mixing of the expected types. In 
the simplest explanation, the deficiency is due to the oscillation of v„ into u T , with Am? ~ 10 -3 eV 
2 and large angle vaccum mixing (sin 2 29 > 0.8) between the two neutrino species. Such oscillation 
is also consistent with the results from the SOUDAN Q and MACRO || experiments. Side by 
side, one also faces the requirement of accounting for the solar neutrino deficit [||. Using the still 
available explanations in terms of v e — oscillation, this means a mass-squared difference of about 
10~ 5 — 1CU 6 eV 2 for the MSW solution (with either small or large angle mixing), or a mass-squared 
splitting smaller by about 4-5 orders for the vacuum oscillation solution. Thus one faces the task 
of reconciling a mass hierarchy with large angle mixing between the two heaviest neutrino species. 
In whichever way that is possible, one has to step outside the domain of the standard electroweak 
model. 

Supersymmetry (SUSY) has been studied for a long time now, from both theoretical and phe- 
nomenological angles, as one of the most attractive options beyond the standard model ||. And 
yet, apart from indirect successes such as offering solutions to the hierarchy puzzle, it has not been 
possible to confront SUSY with any clear experimental results so far. It is therefore appropriate 
that when the atmospheric neutrino results are so emphatically underlining the existence of neu- 
trino masses and mixing with a given pattern, the relevance of SUSY in generating and explaining 
such a pattern should be thoroughly explored. 

If the SUSY extension of the standard model (SM) with the minimal particle content has to 
be invoked for the purpose, then the absence of any right-handed neutrino (and therefore of Dirac 
masses) in the scenario immediately points towards Majorana neutrinos as the likely solution. 
However, the latter implies the violation of lepton number by 2 units. That again can come rather 
naturally in a scenario where R-parity, defined as R = (_i)3B+£+2S^ j g violated. 0. In such a 
scenario, the AC = 1 terms in the Lagrangian can ultimately give rise to Majorana masses either 
through a tree-level see-saw type mechanism or via loop effects. The 'attractive' point here is 
that one does not have to postulate the existence of any particle specifically for the generation of 
neutrino masses, since the superparticles that are indispensable components of a SUSY theory are 
sufficient for the purpose. 

Perhaps the most convenient and cogent (though not unique) way of introducing R-parity 
violation is an extension of the superpotential, using trinilear and/or bilinear terms. Trilinear terms 
in the superpotential Q give rise to L (or i3)-violating Yukawa-type interactions and trilinear soft 
terms in the scalar potential. The presence of bilinear /^-violating terms ||, [9|, [l(J, on the other 
hand, are generally responsible for non-vanishing vacuum expectation values (vev) for sneutrinos, 
which also lead to mixing between neutrinos and neutralinos as also between charged leptons 
and charginos (and similarly between the Higgs and charged slepton/sneutrino states in the scalar 



sector). The bearing of both approaches on neutrino masses have been studied extensively in recent 
times, whereby constraints on the R-parity violating parameters from neutrino masses have been 
derived [11, [ll|]. One also finds in the literature discussions on how to test the consequences of the 



corresponding scenarios in accelerator-based experiments |[3|]. 

In several of these references, it was shown how one could accommodate the mass hierarchy 
together with large angle mixing between and v T through neutrino-neutralino mixing via the 
bilinear terms mentioned above, while the smaller mass splitting between the and v e could 
be due to loop-induced effects. However, one can still ask a number of questions related to the 
parameter space of the theory, extending both over the fermionic and the scalar sector, before the 
SUSY explanation of the SK results can acquire sufficient credibility. In this paper we have tried 
to find answers to some such questions, and to establish that the solution space for the SK deficits 
is not a fine-tuned one. 

To be more specific, our analysis includes both the scalar and spin-1/2 sectors of an R-parity 
violating scenario (where the simplest, bilinear R- violating terms are introduced as necessary in- 
gredients but no generality is otherwise discarded) which can explain the SK data. Keeping the 
value of AmJ T and the large mixing angle as inputs, we have tried to find as general answers as 
possible to the following questions: 

• Is it a necessity to have a large hierarchy between the bilinear terms corresponding to different 
lepton families in the superpotential? 

• Is it required to have a hierarchy between soft /^-violating terms involving different lepton 
flavours? 

• How crucial is it to have a hierarchy of orders between the /^-conserving and /^-violating 
bilinear terms in the superpotential and the scalar potential? 

• How is the suppression of flavour changing neutral currents (FCNC) ensured in such a picture? 

In this analysis, we have used low-energy values of all the parameters in the theory, and have 
not attempted to link them to any specific high-scale physics. However, our chosen convention 
of writing the soft bilinear terms in the scalar potential carries some influence of a supergravity 
(SUGRA)-based model where such terms can be shown to originate from an interference of terms 
belonging to the hidden and observable sectors in the superpotential. 

Section 2 sets up the general framework within which we operate, incorporating the fermion and 
scalar mixing schemes. The different choices of basis in which we have worked are specified there, 
and the parameters which we ultimately use to explore large angle mixing in the most general sense 
are defined. Section 3 contains numerical studies of the constraints on the theory in terms of these 
parameters. This gives us full or partial answers to the questions listed above, obtained either from 
detailed calculations or from simple estimates of order. We conclude in section 4. 



2 The Framework 



The MSSM superpotential is given by (suppressing the SU{2) indices) 

Wmssm = ^HiH 2 + h\ 3 UH x E] + h%QiH x D] + h%QiH 2 U$ (1) 

where /i is the Higgsino mass parameter and the last three terms give all the Yukawa interactions. 
When R-parity is violated, the following additional terms can be added to the superpotential: 

Wft = A,, 7, /.,/., /--X + \jkLiQjD% + Xiji.l 'i + EiLiH.2 (2) 

Where the A"-terms correspond to ^-violation, and the remaining ones, to ^-violation. The 
absence of proton decay makes it customary to have one of these two types of nonconservation at 
a time. In the rest of this discussion, we shall not consider S-violation. 

The A-and A'-type terms have been widely studied; their contributions to neutrino masses can 
be only through loops, and their multitude (there are 36 of them altogether) makes the necessary 
adjustments possible for creating the required values of neutrino masses and mixing angles. 

More interesting, however, are the three bilinear terms eiL{H 2 . It is in fact quite useful to use 
them as the starting inputs of R-parity violation in the theory. There being at most only three 
parameters of this type, the model looks much more predictive than one with 36 unrelated trilinear 
terms. Furthermore, the physical effects of the trilinear terms can be generated from the bilinears, 
by going to the appropriate bases. There are additional interesting consequences of the bilinear 
terms. The presence of the Li^-term means a mixing between the Higgsinos and charged and 
neutral lepton states. In addition, the scalar potential in such a case contains terms bilinear in the 
slepton and the higgs fields and involving only second and third generation of sleptons (the reason 
behind this assumption will be clarified as we proceed) the terms are as follows : 



V sca i = m 2 Ls L 2 3 + m\J\ + m 2 H 2 + m 2 2 H 2 + B^H X H 2 
+ B 2 e 2 L 2 H 2 + B 3 e 3 L 3 H 2 + [ie 3 L 3 B\ + /j,e 2 L 2 Hi + (3) 

where denotes the mass of the i-th scalar doublet, Li = I _ I , i being 2 and 3 for the second 

and third generations respectively, at the electroweak scale. Here mx 2 and rriL 3 are the slepton 
mass parameters. In our subsequent analysis, all the sleptons (of both chiralities) and sneutrinos 
have been assumed to be degenerate. 

An immediate consequence of the additional (£- violating) soft terms in the scalar potential is 



a set of non- vanishing sneutrino vev's [14|. This is a characteristic feature of this scenario, which 
in addition produces neutrino (charged lepton)-gaugino mixing via the sneutrino-neutrino(charged 
lepton)-gaugino interaction terms. The former type of mixing leads to a neutrino mass at the 
tree-level. 



Since our primary goal is to explain large angle — v T mixing, we simplify the picture by 
assuming that only the second and third generations enter into this tree-level mixing process. For 
this, we postulate only two bilinear R-violating terms, proportional to e 2 and €3. These two terms, 
together with the soft bilinear terms B 2 and B 3 , form the set of independent R-parity violating 
parameters in this basis, henceforth to be called basis 1. The vev's corresponding to the muonic 
and tau sneutrinos are and v T respectively in this basis. For reasons that will become apparent 
later, we choose to treat these vev's as independent parameters, and use them to derive values of 
the soft terms. For that purpose, one has to make use of the set of tadpole equations arising from 
electroweak symmetry breaking: 



(ml + 2\c)vi + Bfiv 2 + /xe 2 f M + fie 3 v T = 
(m\ - 2Xc)v 2 + B/ivi + B 2 e 2 v pi + B 3 e 3 v T = 
(m% + 2Ac)f M + B 2 e 2 v 2 + fie 2 v 1 + e 2 e 3 v T = 
(m\ + 2Xc)v T + B 3 e 3 v 2 + fie 3 v 1 + e 2 e 3 v^ = 



(4) 
(5) 
(6) 
(7) 



where v\ =< H\ >, v 2 =< H 2 >, c 



)\ - v\ + v 2 T + vl) and A = (g 2 + g' 2 )/8. 



While two of these equations can be used to eliminate the soft Higgs mass terms m\ and m 2 , 
the /^-violating soft terms B 2 and B 3 can be obtained from the remaining two: 

1 



Bo 



B, = — 



£2^2 
1 

£3^2 



(e 2 e 3 v T + 2Xcv^ + m^v^ + e 2 ^,v 1 ) 
(e2e 3 Vfj, + 2Xcv T + m\ T v T + e 3 fivi) 



(8) 
(9) 



where, again, the sneutrino masses have been assumed to be degenerate with a common slepton 
mass parameter. 

The next step is to rotate away both the e-terms from the superpotential. In the process we go 
from the basis (Hi, L 3 , L 2 ) to (H[, L f 3 , L' 2 ) using the following rotation : 



L' 3 
\L' 2 J 



( 



where s 2 



V4+> 



c 2 



£3 



2' 



= $ 



c 3 s 3 c 2 S 3 S 2 
-s 3 C 3 C 2 C 3 S 2 

-s 2 c 2 J 

V4+4 



\ ( Hi\ 

L 3 
L 2 



(10) 



V ^ 2 / 



«3 



, and // = J fi 2 + e\ + ei. Clearly, 



this leaves ii!H' x H 2 as the only bilinear term in the superpotential. The physical consequences 
of bilinears R-parity violation, however, are still existent in this basis, since the scalar potential 
contains its signature, and one has now 'rotated' sneutrino vev's v'^ and v' T which, together with 
the Higgs vev v[, are connected to the set (vi,v T ,Vn) by the rotation matrix given in Eqn. (10) 
above. 

In this basis (called basis 2), the neutralino mass matrix (which is of dimension 6x6 after 
including the neutrinos) is given by 
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where the successive rows and columns correspond to (H2, H[, —iW 3 , —iB, v' T , 1/ ). Here M and M' 
are the SU(2) and U(l) gaugino mass parameters respectively, // is the Higgsino mass parameter 
in basis 2 and v[ =< H[ >. This leads to one non- vanishing tree-level neutrino mass eigenvalue. 
The root of this can be traced to the fact that one linear combination of u' and u' T , given by 



1/3 = u T cos 6 + v ' sin 9 (12) 
enters into cross-terms with the B and W3, while its orthogonal state V2 decouples from the mass 



matrix and therefore remains massless 15]. The angle 6 is given by 



cos 9 



v'Jv' 



(13) 



with v' 



+ v' T . The quantity v' , which is a kind of 'effective' sneutrino vev in a basis 



where the bilinear terms are rotated away from the superpotential, is a basis-independent measure 
of R-parity violation, which directly controls the tree-level neutrino mass acquired in the process. 
The mass implied by the atmospheric i/^-deficit in the SK data requires v' to be in the range 
(1 — 3) x 10~ 4 GeV approximately |l(|, depending on the mass of the lightest neutralino (x?)- It 
is also interesting to note at this stage that the smallness of the neutrino mass does not limit the 
value of the e-parameters so long as v' lies within about 100 keV or so. It has been demonstrated, 
for example, that in theories based on N=l supergravity (SUGRA), it is indeed possible to have a 
small v' in spite of large values of the e's, by setting all bilinear soft terms (both /^-violating and 



conserving) to the same value at the SUGRA breaking scale [17]. 

Modulo the very small neutrino-neutralino mixing, the angle 9 defined above can be identified 
with the neutrino vacuum mixing angle (0q) if the charged lepton mass matrix is diagonal in basis 
2. Only in such a case can we say that a near-equality of the vev's v' and v' T is required by the 
condition of large-angle mixing. However, the general form of the charged lepton mass matrix in 
this basis is 



Mi 



(14) 



/3(-^S 2 S 3 + C 2 V' 1 ) f2{s 2 v' 1 + C 2 S 3 v' fl ) 
f3(-C3S 2 v' 1 + S 2 S 3 v' T ) f2(c2C 3 v[ - C 2 S 3 v' T ) 

Here we have assigned (t' l , along the rows and (tr, hr) along the columns. One should note 



that h = h 33 = ^ and f 2 = h\ 



'22 



. Thus Mi has non- vanishing off-diagonal terms in general. 



The conditions under which it can be approximately treated as diagonal (and the equality of the 
two sneutrino vev's in basis 2 is a necessity) will be specified in the next section. Instead, we 
observe here that the actual neutrino mixing matrix is given by 

Vi = V e U T (15) 

where Vq is the 2x2 orthogonal matrix corresponding to the rotation angle 9, and U is the matrix 
that diagonalises Mi is basis 2. It is the matrix Vi which should correspond to a rotation angle of 
7r/4 for maximal mixing, and to a range approximately between 32 and 58 degrees for sin 2 26q > 0.8, 
#o being the neutrino mixing angle determined by V\ Q. This restricts one to specific allowed ranges 
in the ratio of the vev's v'^ and v' T , for each value of the ratio €2/63. Corresponding to these allowed 
values, the 2-dimensional space in the soft parameters B2 and B3 also gets constrained. In the next 
section, we shall present a detailed map of the region of the parameter space which corresponds to 
large angle mixing as depicted by the SK data. 

Before we end this section, we also list the three scalar mass matrices for this theory in basis 
1. These include both the Higgs and the slepton/sneutrino sectors, with the appropriate kinds of 
mixing between them. 

The charged scalar mass matrix is given by 
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(16) 



with 
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1 This is the case if the angle is to lie in the first quadrant. As we shall see in detail in the next section, there can 
be nontrivial solutions for 8q in the second quadrant also. 
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I 2,2 2 , 2\ 
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where both the left-arid the right-chiral sleptons for each flavour have been included, the basis being 
(Hi, H2, tl, jxi,, tr, flu). Similarly, in the bases (Re(Hi), Re(H2), Re(v T ), Re(Vfj)) and (Im(H{), 
Im(H2), Im(u T ), Im(ua)) respectively, the neutral scalar and the pseudoscalar mass matrices are 
given by 
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(17) 



and 



Mi 
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-B 2 e 2 
e 2 e 3 
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(18) 



/ m\ + 2Ac — By ye 3 ye2 \ 

—By, m\ — 2Ac —-6363 

ye 3 -B 3 e 3 m 2 ~ T + 2Ac 

ye 2 -B 2 e 2 e 2 e 3 m?^ + 2Ac / 

It may be remarked that the charged scalar and the neutral pseudoscalar mass matrices will 
each have a zero eigenvalue, corresponding to the Goldstone bosons. Also, these scalar matrices 
have to be used for a complete determination of the allowed space in B 2 and B 3 , since a number 
of conditions related to electroweak symmetry breaking need to be satisfied before values of these 
parameters, as extracted from Eqns. (8) and (9), pass off as valid ones. These conditions |18|1 
include the requirement that the potential be bounded from below, the negativity of the vev's and 
non- negativity of the eigenvalues. They can be obtained by a straightforward generalisation of the 
corresponding conditions given in reference 18 where lepton number violation in only one family 
has been considered. 



3 The constraints 



As we have already mentioned, our first constraint is on the quantity v' . The allowed range of 
Am^ T , combining the fully contained events, partially contained events and upward-going muons, 



is about 1.5 — 6.0 x 10 3 eV 2 at 90% confidence level [19|. For the lightest neutralino mass varying 
between 50 and 200 GeV, this corresponds to v = 0.0001 - 0.0003 GeV (100-300 keV) to a fair 
degree of approximation. Since v = \Jv'^ + v'^, this automatically puts a constraint on v'^ and v' T , 
and hence on the vev's in basis 1 for each value of the e-parameters. 

Next, the value of the angle 9 should determine v'^ and v' T completely, once v is fixed. As we 
have pointed out in the previous section, for each value of £2/63, one is restricted to particular 
values of 6 in order that the condition of large angle mixing, namely sin 2 29q > 0.8, is satified in 
terms of the ultimate neutrino mixing angle 9q. In Figs. 1 and 2 we outline these allowed regions 
in the parameter space of €2/63 vs. 9, for €2/^3 < 1 and €2/63 > 1 respectively. The allowed bands 
are found to be sensitive to the ratio rather than the actual values of the e's. The interesting point 
to note here is that there are two values of 9 for each value of the ratio along the x-axis. This 
is due to the two solutions with angles in either of the first or the second quadrant, both cases 
yielding the same value for the oscillation probability. Physically, the second case corresponds to 
the superposition of the and v T states being performed with an extra phase rotation through 
an angle ir for one of them. In other words, the two solutions represent situations with the two 
neutrino flavours having identical and opposite CP-properties respectively. 

The angles 9 and #0 are practically equal if || << In this case the off-diagonal elements of 
the charged lepton mass matrix Mi in basis 2 are neligibly small compared to the diagonal ones. 
However, Figs. 1 and 2 clearly demonstrate that such a large hierarchy between the two bilinear 
terms in the superpotential is by no means a necessity and that in fact it is possible to have them 
not only with the same order of magnitude but also with an inverted hierachy as well. 

Each point in the allowed regions shown in Figs. 1 and 2 corresponds to a theory for which the 
two bilinear soft terms B2 and B3 can be calculated. The ranges of values thus obtained have been 
plotted in Figs. 3 and 4, again with opposite hierarchies of the parameters £2 and €3. The two 
bands in each figure arise from the two solutions of 9 for each ratio of the e's and a particular value 
of €3 (Fig. 1) or €2 (Fig. 2). This underlines the fact that the two solutions for the angle in the 
two quadrants are indeed physically distinct, with all the different phenomenological implications 
of the two combinations of the £-violating soft terms 0. 

The value of B2 in Fig. 3 is found to become very large when €2 is assigned a value much smaller 
than that of €3 (and similarly for B3 in Fig. 4 with a very small €3). This is analogous to the 
problem of the MSSM. However, Figs. 1 and 2 already tell us that such a large hierarchy of the 
2 In fact, there could be two more sets of solutions in the two remaining quadrants. However, these solutions are 
not distinct from the previous ones, as they can be mapped back to the cases of same and opposite parity for the two 
neutrinos. Explicit computaions also lead to the two already obtained combinations of B2 and B3. 



two e's is not a necessary condition for large angle neutrino mixing. The very expressions for the 
.B-parameters show that except in cases where very small e 2 ( 3 ) jacks them up to large values, they 
are also controlled by the scale of the slepton mass. We also see from Fig. 3 that B3 cannot attain 
values much larger than the electroweak scale so long as e 2 < £3. The same feature is observed for 
B 2 in Fig. 4. 

However, it must be noted at this point that the sets of values for £> 2 and B% as seen in Fig. 3 
or 4 result from the choices of two parameters, viz., e 3 (Fig.3) or e 2 (Fig.4) and their appropriate 
ratios, the choices being rather specific and spanned over identical regions in these two figures. 
This is the reason why one of B2 and B% ranges over several orders of magnitude while the other 
is highly constrained in a particular figure. Here, it should be made clear that the magnitudes of 
B2 and B3 are mainly determined by those of e 2 and €3 respectively for very small values (< 1CT 5 ) 
of the latter two. As, neither the smallness of the individual values of e's nor that of their ratios is 
anyway restricted, we can always tune these two parameters so that both of B2 (Eqn. (8)) and B3 
(Eqn. (9)) are simultaneously very large (compared to the electroweak scale) leading to additional 
regions in the allowed 5 2 — B3 space which are absent in Figs. 3 and 4. 

Summing all these up, the parameter space of the soft terms (modulo their additional depen- 
dence on the MSSM parameters) is seen to be sufficiently free from any requirement of fine-tuning 
or large hierarchy, and large angle neutrino mixing as well as the expected neutrino mass hierarchy 
in the second and third families is reproduced over a rather wide range in the space of R-parity 
violating parameters. Also, B2 and B3 are not compelled to show any hierarchical behaviours with 
respect to the R-conserving soft term B, and any approach connecting them at a high scale is 
consistent with the range of values allowed here. 

The discussion becomes more transparent if we break up the results shown in Figs. 3 and 4 into 
3 broad regions, viz., (i) very small e 2 (e3) (~ 10~ 8 ), (ii) somewhat intermediate values of e 2 (es) 
(~ 10~ 5 ) and (iii) rather large values of 62(63) (~ 1CT 2 ) in reference to Eqn. (8) (Eqn. (9)). The 
numbers in the parenthesis correspond to the left extreme graphs of Fig. 3 (Fig. 4). In case (i) 
£2(^3) dependent terms lead and its presence in the denominators with a very small value leads 
to a large value for i? 2 (i?3). When e 2 (e3) increases as for case (ii) the contributions from these 
terms gradually become comparable to that coming from the 4th term which depends only upon 
MSSM parametrs, viz., ll and tan/3. This restricts i? 2 (i?3) to intermediate values. The relative 
sign between the 4th term and other terms collectively is instrumental in fixing the sign of i? 2 (£?3) 
in this region and is clearly visible as we proceed along the second row of graphs in Fig. 3 (Fig. 
4). In case (iii) the 4-th term controls B 2 (Bs). Naturally, in this case, the order of ^2(^3) is set 
by \i and tan (3 and is restricted to be around 100 GeV for choices of these two parameters that are 
compatible with existing collider data. 

In the results presented above, fx and ll have been taken to be of the same sign, with both e 2 
and €3 having signs opposite to it. On reversing this relative sign (and also that between the two 
e-terms), it is seen that the orders of magnitude of B2 and B3 do not change; however, the signs of 



one or both of them are liable to get reversed. 

On reversing the sign of /u, both B2 and B3 may change sign if the bulk contributions to them 
come from the 4-th terms in Eqns. (8) and (9). Otherwise, the relative sign between B2 and B% 
and also their magnitudes (to a good approximation) are preserved under such reversal. Changing 
the sign of \x should be accompanied by a change in sign of B (of MSSM) to retain a relative sign 
between themselves which is a must to achieve electroweak symmetry breaking at the proper scale. 
On the other hand, B2 and B3 pick up a relative sign on allowing for the same between €2 and €3 
only if the latter two contribute heavily to B2 and B3 respectively. 

We restrict ourselves by presenting only two sets of figures (Figs. 3 and 4) which jointly cover 
rather wide intervals of €2 and €3. Results of any combinatorics of relative signs between \i and the 
e's can be directly estimated to a fair degree of accuracy from these two sets of figures. These show 
that the magnitudes of B2 and B3 can range over several orders of magnitude with any relative 
sign between them. 

It should be re-iterated that variations of B2 (-B3) with the MSSM parameters fi and tan/3 are 
likely to be most pronounced when major contributions to it comes from the 4th term in Eqn. (8) 
(Eqn. (9)). This should be contrasted with the results presented in Figs. 1 and 2 where the allowed 
band is quite insensitive to the MSSM parameters. B2 and B3 exhibits such a sensitivity when €2 
(€3) is ~ 10~ 3 or more, so that the other terms in Eqn. (8) (Eqn. (9)) have comparable orders 
of magnitude (~ electroweak scale) leading to B2 (-B3) lying around 100 GeV. Here, variation of 
tan/3 over a range 5 to 50 can change B2 (-B3) at most by 1 order while variation with \i over a 
range —500 to +500 does not change things much, although a flip in the relative sign between B2 
and -B3 may be a possible outcome of the variation of [i. 

In the above, we have not discussed the mass-splitting required between the electron and the 
muon neutrinos for explaining the solar neutrino problem. As we have already mentioned, the 
trilinear interactions in Eqn. (2) give rise to loop-induced masses that can account for this splitting 
rather well. The nature of the splitting will depend on whether one wants to provide an explanation 
in terms of matter-enhanced or vacuum oscillations, and the only consequent constraints will come 
on the respective A-and A -couplings in basis 2. 

There are, however, additional issues that must be addressed here. In general, the very structure 
of the charged scalar, neutral scalar and pseudoscalar mass terms, obtainable from Eqn. (3), admits 
of mixing between the smuon and stau flavour states, as also between the corresponding sneutrino 
flavours in both the neutral scalar and pseudoscalar mass matrices. This in general causes a 
mismatch between these scalar mass matrices and the charged lepton mass matrix, standing in the 
way of their simultaneous diagonalisation. This in general can lead to an enhancement of FCNC, 
contributing to processes like r — > ^7. 

FCNC suppression can in general be ensured by one of three mechanisms- degeneracy, alignment 
and decoupling pQ| . Decoupling in this case requires the slepton and sneutrino masses to be very 
high so that contributions to the relevant processes from loop diagrams can be suppressed. With 



these masses on the order of 100 GeV, one has to ensure either alignment of the scalar mass matrices 
with the one for charged leptons, or a close degeneracy of the scalar eigenvalues. Both conditions 

are seen to be satisfied in this case if (M~ - << (M~ ~ )i for i = s,p or c (i.e. in the neutral 

L2-L3 L2L2 

scalar, pseudoscalar or charged scalar mass squared matrices). 

Figures 5 and 6 contain plots of the ratio of the ultl and /j,l/J-l elements for each of the three 
matrices, against the ratio of £2 and €3 (again for both £2 < £3 and £2 > £3). A pattern is revealed 
in each case on imposing the requirement on the off-diagonal elements to be small compared to the 
diagonal ones. 

First consider the case £2 < £3- For small values of this ratio, the smallness of all the diagonal 
terms (and therefore the suppression of FCNC) is assured, irrespective of the value of £3. In such 
cases it is enough to consider the scalar mass matrices in basis 2 itself, since the condition £2 < < £3 
makes the charged lepton mass matrix practically diagonal in this basis. The off-diagonal terms 
in this basis are controlled by the quantity v' 2 which is constrained from neutrino masses to be 
much smaller than m~ and /U 2 which dominate the values of the diagonal terms. However, as one 
increases £2 A3, the charged lepton mass matrix ceases to be diagonal in basis 2, and the scalar 
mass matrices need to be evaluated after diagonalising it. In such cases the off-diagonal elements 
are also influenced by the actual value of £3. As can be seen from Fig. 5, for £3 approaching 100 
GeV, the off-diagonal term starts becoming comparable to the diagonal ones as the ratio £2/f3 is 
close to unity. Exactly the same thing can be said of £2 in Fig. 6 where predictions for an inverted 
ratio between the two e's are shown. 

A rather interesting conclusion follows from the above discussion. While a compatibility with 
the SK results is generally observed over a large area in the space of R-parity violating parameters, 
the need to suppress FCNC effects introduces at least one type of hierarchy. This is apparent 
from the fact that one of the two bilinear parameters £2 and £3 can be allowed to be as large as in 
the electroweak scale if their mutual ratio has a hierarchy of 2 to 3 orders of magnitude. On the 
other hand, allowing the two of them to be of the same order is consistent with the suppression of 
FCNC only if the scale of their magnitude is small compared with the corresponding R-conserving 
parameter i.e. £2(3) « [i. Thus, cither a hierarchy in the two R- violating parameters in the 
superpotential or a smallness of the R- violating ones compared to R-conserving ones is a necessary 
consequence of this scenario. 

4 Conclusion 

We have investigated a SUSY theory where the bilinear R-parity violating terms are responsible 
for tree-level neutrino mass and large angle mixing, leading to a pattern that explains the observed 
defhciency of atmospheric muonic neutrinos. The scenario also admits of trilinear R-parity violating 
terms that cause the mass splitting necesary to explain the solar neutrino puzzle. 

A detailed study of the parameter space of the theory shows that, as far as constraints arising 



from the mass-splitting and large angle mixing are concerned, neither a hierarchy of R- violating 
parameters nor a fine-tuning among them is necassary, and a large area of the parameter space of 
the R-breaking soft terms is 'naturally' allowed. However, the suppression of FCNC can require 
either the R-violating parameters in the superpotential to be small compared to the //-parameter, 
or such parameters themselves to have a mutual hierarchy of values. 
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Figure 1: Allowed bands of neutrino mixing angle (0 in degrees) corresponding to maximal mixing 
between 2nd and 3rd generation of neutrinos are plotted against the ratio of epsilons (f|). The 
MSSM parameters are \i = —500, tan (3 = 10 and B = 100. The pattern remains almost the same 
for different es's. The different shadings arise due to varying number of sample points in different 
regions. 
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Figure 2: Allowed bands of neutrino mixing angle (6 in degrees) corresponding to maximal mixing 
between 2nd and 3rd generation of neutrinos is plotted against the ratio of epsilons (||). The 
MSSM parameters are same as in Fig. 1. The pattern remains almost the same for different e2's. 
The different shadings arise due to varying number of sample points in different regions. 
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Figure 3: Allowed regions of Bi - B% space (in GeV). R = f| varies along the columns as indicated 
in the top margin. £3 varies along the rows as shown in the right margin. The MSSM parameters 
are n = —500, tan/3 = 10 and the common slepton mass at the electroweak scale is 200 GeV. 
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Figure 4: Allowed regions of B2 - B% space (in GeV). Ri = || varies along the columns as indicated 
in the top margin. 62 varies along the rows as shown in the right margin. The MSSM parameters 
are same as in Fig. 3. 



M 2 [3,4] 
M 2 [4,4] 




0.0001 0.001 0.01 0.1 1 

Figure 5: Ratio of /Ultl ([3,4]) and hl^l ([4,4]) terms for different mass-squared matrices (in a basis 
where the lepton mass matrix is diagonal) is plotted against the ratio (f|). The MSSM parameters 
are n = —500, tan /? = 10, B = 100 and common slepton mass at the electroweak scale is 200 GeV. 



le-15 



1 1 — 1 1 1 1 II 1 1 1 — 1 1 1 1 1 1 1 1 1 — 1 1 1 1 1 1 1 1 1 — 1 1 1 1 1 1 

-t r\r\ r— — ~-*~ =r ^* r ~^^ era 






pseudo-scalar 

cnarged scalar _ 






e 2 = 0.0001 




i i 





0.0001 0.001 0.01 0.1 1 

f-1 

Figure 6: Ratio of \xltl ([3,4]) and hl^l ([4,4]) terms for different mass-squared matrices (in a basis 
where the lepton mass matrix is diagonal) is plotted against the ratio (^). The MSSM parameters 



are same as in Fig. 5. 



